A Numerical Scheme for the Pore-Scale Simulation of Crystal Dissolution and Precipitation in Porous Media

In this paper we analyze a numerical scheme for a dissolution and precipitation in porous media. We focus here on the chemistry, which is modeled by a parabolic problem that is coupled through the boundary conditions to an ordinary differential inclusion defined on the boundary. We use a regularization approach for constructing a semi-implicit scheme that is stable and convergent. For dealing with the emerging time discrete nonlinear problems, we propose a simple fixed point iterative procedure. The paper is concluded by numerical results.